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^ : Abstract. Basing on the increasing evidence for the cosmological relevance of the local Hubble flow, we consider a 

simple analytical cosmological model for the Local Universe. This is a non-Friedmann model with a non-uniform 
static space-time. The major dynamical factor controlling the local expansion is the antigravity produced by the 
omnipresent and permanent dark energy of the cosmic vacuum (or the cosmological constant). The antigravity 
dominates at distances larger than 1-2 Mpc from the center of the Local group. The model gives a natural 
, explanation of the two key quantitative characteristics of the local expansion flow, which are the local Hubble 

constant and the velocity dispersion of the flow. The observed kinematical similarity of the local and global flows of 
expansion is clarified by the model. We demonstrate analytically the efficiency of the vacuum cooling mechanism 
that allows one to see the Hubble law so close to the Local group. Special significance is argued for the 'universal 
^ ' Hubble constant' Hy, depending only on the vacuum density {Hy — 60 km/s/Mpc). The model makes a number 

i—{ I of testable predictions. It also shows unexpectedly that it is the dwarf galaxies of the local flow with the shortest 

. distances and lowest redshifts that may be the most sensitive indicators of dark energy in our neighborhood. 

Key words, dark matter - cosmological parameters - Local Group 

1. Introduction so much that, for the first time, we have reliable grounds 

for finding a solution to the Hubble-Sandage paradox. The 
Cosmology has changed dramatically with the discovery paradox is in the fact that cosmological expansion with 
of the cosmic acceleration and dark energy of the cosmic a regular linear velocity- distance relation was discovered 
vacuum (or the cosmological constant) from observations by Edwin Hubble where it should not be. Indeed, the 
of high-redshift supernovae (Riess et al. 1998; Perlmutter Friedmann theory predicts and describes the expansion of 
et al. 1999). A new 'concordance cosmological model' has a smooth, homogeneous self-gravitating medium. But lo- 
emerged which describes the accelerating cosmological ex- cally, the matter is far from uniform and rather reveals a 
pansion on cosmological scales beyond the cell of unifor- fairly chaotic spatial distribution of galaxies. Therefore the 
mity (> 300 - 1000 Mpc). The model is in good agree- Friedmann theory cannot be apphed to the Local Volume, 
ment with all cosmological data including the WMAP re- The paradox was splendidly realized by Allan Sandage 
suits on the CMB anisotropy (Spergel et al. 2003). The (Sandage et al. 1972, Sandage 1986, 1999). He saw a most 
Friedmann theory with a non-zero positive cosmological striking mystery in the fact that the local expansion pro- 
constant serves as a reliable basis for the concordance ceeds in a quite regular way: the linear velocity-distance 
model, the present-day standard model in cosmology (see relation (the Hubble law) with a low velocity dispersion is 
e.g. reviews by Chernin 2001; Peebles & Ratra 2003). clearly observed here. Moreover, the local rate of expan- 
A new field of studies in which the cosmic vacuum has ^io" is similar, if not identical, to the global rate charac- 
proved to be important, lies deep inside the cell of uni- Prized by the Hubble constant. 

formity. In this Local Volume (distances < 10 — 20 Mpc) Recently, Thim et al. (2003) have emphasized a related 

cosmological expansion was originally discovered. Our un- problem: there is no evidence for the expected turn-down 

derstanding of the physical conditions there has improved of the velocity-distance relation, when one approaches the 
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Local Group (c.f. Ekholm et al. 2001). Whiting (2005) 
notes the similar absence of infall flow around local galaxy 
groups in general. 

Soon after the discovery of the cosmic vacuum, we sug- 
gested (Chernin, Teerikorpi, Baryshev 2000 - hereafter 
Paper I) that the cosmic vacuum with its omnipresent and 
perfectly uniform density might provide the dynamical 
background for a regular quiescent expansion in the Local 
Volume. Our key argument comes from the fact that the 
vacuum antigravity dominates over the matter gravity of 
the Local Group at the distances larger than 1-2 Mpc from 
the center of the group (Chernin 2001, Baryshev, Chernin, 
Teerikorpi 2001, Karachentsev, Chernin, Teerikorpi 2003). 
The increasing observational evidence and theoretical ar- 
guments for this view have been recently discussed by 
Teerikorpi, Chernin, Baryshev (2005). In that paper it 
was also shown that dark energy may explain still an- 
other problem in the Local Volume: why peculiar veloci- 
ties do not depend on luminosity (mass) (Karachentsev & 
Makarov 1996; Whiting 2003; 2005). 

This new approach has extended cosmology into the 
deep interior of the cosmic cell of uniformity. Especially, 
the Local Volume must be regarded as the nearest cos- 
mologically relevant object. So we find it good to use the 
term 'the Local Universe' for our relatively small 'Local 
Volume' sample of the Universe which, however, is fairly 
typical for galaxy space in general. We will also speak 
about 'local cosmology', assuming at the same time that 
the Universe as a whole is described by the concordance 
model of the traditional 'global' cosmology. 

In the present paper, we develop our approach fur- 
ther and study a simple analytical non-Priedmann model 
for the Local Universe (Sees. 2, 3). The model naturally 
explains two key characteristics of the local expansion 
flow: the local Hubble constant and the velocity disper- 
sion (Sec.4). A special significance is argued for the uni- 
versal 'vacuum' Hubble constant (Sec. 5). The model pre- 
dicts that the cosmic antigravity and vacuum density can 
be detected and measured in our nearest neighbourhood 
(Sec. 6). In Sec. 7, complementary computer simulations 
are briefly reviewed and a general discussion of the results 
is given. We summarize the conclusions are in Sec. 8. 

2. Local cosmology: Newtonian theory 

It is clear from the considerations of Sec.l that a 
non-Priedmann counter-part of Priedmann's cosmological 
model is needed for the Local Universe where the mat- 
ter distribution is highly non-uniform. In constructing a 
non-Priedmann theory, we focus first on the nearest and 
most important distances, 1 — 5 Mpc, where the obser- 
vational data are fairly complete and distance determi- 
nations are accurate enough (Karachentsev & Makarov 
2001, Karachentsev et al. 2002, 2003a, Thim et al. 2003, 
Karachentsev et al. 2004; Paturel & Teerikorpi 2005, 
Rekola et al. 2005). We start with a simple version of 
the theory in which the following physical conditions are 
assumed for the present epoch of cosmic evolution: 



1. The environment of the Local Universe (i.e. the dis- 
tribution and motion of matter outside it) does not affect 
essentially the dynamics inside it and one may neglect the 
matter distribution and motions at the distances larger 
than, say, 5 Mpc, in the first and main approximation. 

2. The matter mass of the Local Universe is collected 
almost entirely - with a few percent accuracy - within the 
area of the Local Group at the distances < 1 Mpc from 
the barycenter of the group, in this approximation. 

3. The gravity field produced by the Local Group is 
spherically symmetrical and static at the distances > 1.5 — 
2 Mpc from the group barycenter (more about this below). 

4. The interaction of the out-flowing dwarf galaxies 
with each other is negligible. 

5. The flow consists of the galaxies which have escaped 
from the gravitational potential well of the Local Group. 

2.1. The equation of motion 

Under these simplified, but for the present discussion suffi- 
ciently reasonable assumptions, the dynamics at distances 
1-5 Mpc from us is reduced to the problem of a test parti- 
cle motion in the gravity field of the Local Group treated 
as a point-like mass on the background of the antigrav- 
iting vacuum. Gravity of the Local Group produces the 
radial force — GM/r^, where r is the distance from the 
group barycenter and M = 1.5 X lO^^M© is the group 
mass. The antigravity of vacuum produces the radial force 
G2pvi^r^)/r^, where — 2py = py + Spy is the effec- 
tive gravitating density of vacuum, according to General 
Relativity (see e.g. Chernin 2001), and pv = 7 x 10~^° 
g/cm'^ is the concordance value of the vacuum density. 

The radial component of motion in this grav- 
ity/antigravity force field obeys the Newtonian equation: 

r{t,x) = -GMy+r/A\ (1) 

where r{t, x) is the distance of a particle to the barycenter 
of the Local group, x is the Lagrangian coordinate of the 
particle; the constant 

A = (^p^)-i/2 ^ 5 X 1017 s 1.5 X lO^^ cm (2) 

o 

is the characteristic vacuum time/length; hereafter the 
speed of light c = 1 is used in the formulas. It is interesting 
that the constant vacuum time/length is numerically near 
the current cosmic age or the distance to the cosmological 
horizon: to ^ A; a, discussion of this fact may be found, 
for instance, in Chernin (2005). 

It is seen from Eq.l that the gravity force (cx 1/r^) 
dominates over the antigravity force (oc r) at smaller dis- 
tances; the acceleration is negative there. At the distance 

Rv = (GMA2)1/3 = (AM/py)i/3 ~ 1.3 Mpc (3) 

the gravity and antigravity balance each other, so the ac- 
celeration is zero at the 'zero-gravity sphere' of the radius 
Rv (Paper I). At larger distances, r > Ry, antigravity 
dominates, and the acceleration is positive. 
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We will study mostly the region outside the zero- 
gravity sphere where the gravity field is assumed to be 
static and spherically symmetrical (see above). In fact, 
computer simulations of the Local group and its gravity 
(Dolgachev et al. 2003, Chernin ct al. 2004) show that 
the zero-gravity surface is fairly spherically symmetrical 
at the present epoch and it is approximately (with a 10- 
15% accuracy, at least) unchanged in size and shape over 
more than 12 Gyr of the history of the Local Group. 

The first integral of Eq.l expresses, as usual, the me- 
chanical energy conservation: 

lr^ = GM/r+\{r/Af+E, (4) 

where E{x) is the total mechanical energy of a particle 
with the Lagrangian coordinate x (per its unit mass). Here 
the potential energy 

U{r) = -GM/r-^{r/Af (5) 

is negative, and both matter and vacuum give to U contri- 
butions of the same sign (contrary to the right-hand side of 
Eq.l). Because of the vacuum, the gravitational potential 
cannot be normalized toU = at i? = oo, and U{R) goes 
to — oo in both limits of i? ^ and R ^ oo which is not a 
standard situation in mechanics. It means, in particular, 
that the trajectories with £^ < are not necessarily finite. 
Such a behaviour of the potential has a clear analogue in 
General Relativity applied to the same problem (sect. 3). 

2.2. Radial trajectories 

The dwarf galaxies of the local expansion flow move along 
approximately radial trajectories in the potential U{r). 
The total energy of each particle that has escaped from 
the gravity potential well of the Local Group exceeds the 
maximal value of the potential U: 

E > Umax = ~^GM/Rv. (6) 

Asymptotically, when r goes to infinity, and the vacuum 
dominates entirely, the solution of Eq.4 is exponential, 

r(i,x) (X expt/A, (7) 

so that the linear velocity-distance law, f = r/Av, appears 
with the constant expansion rate: 

Hv = llA~m km/s/Mpc. (8) 

The quantity Hy is a universal physical constant di- 
rectly related to the vacuum density (or the cosmological 
constant); its provisional numerical value, here adopted, 
comes from the concordance data on py. We call Hy the 
universal Hubble constant. 

Generally, integrating Eq.4, one gets: 

t-T{x)= [ [2GM/r+{r/A)^ + 2E]-^^^dr, (9) 



where T{x) is the second integration constant. Two time- 
independent functions of the Lagrangian coordinates, 
E{x) and fo(x)) are determined by the particle's "initial 
conditions" that are its velocity and distance from the 
center at a given moment of time. 

In the simplest case of the parabolic {E = 0) motion, 
the exact integral of Eq.9 is expressed in terms of elemen- 
tary fimctions: 

r{t,x) oc [sinh{3t/2A)f/^. (10) 
2.3. Local vs. global 

It is significant that the mathematical structure of Eq.4 

is similar to that of the Fricdmann cosmological equation 
for the global Friedmann radius R{t): 

^R{t, xf = GM{x)/R + Ir'/A' + E{x), (11) 

where M(x) is the time-independent mass of the non- 
relativistic matter within the sphere of the radius R{x, t), 
and X is the Lagrangian coordinate which corresponds 
to the Eulerian coordinate R{t,x)', the time-independent 
value E{x) = k{R/a)^, a is the curvature radius, —k = 
— 1, 0, +1 is the sign of the 3D space curvature. As is well- 
known, the function R{t, x) is a product of the curva- 
ture radius a(t) and a function of the Lagrangian coor- 
dinate only. In terms of physics, the difference of Eq.ll 
from Eq.4 is in the dependence of the mass M(x) on the 
Lagrangian coordinate in the Friedmann equation, due to 
the uniform matter distribution. According to the non- 
Priedmann Eq.4, each body moves in the potential of the 
same gravitational mass. Clearly, the solution of Eq.ll has 
the same - in mathematical structure - form as Eq.lO for 
the parabolic {E = 0) motion. 

It is especially suggestive that, in the limit of large R 
and large r, both equations lead to the same exponen- 
tial solution with the linear velocity-distance relation and 
the universal Hubble constant Hy — 1/A. This similarity 
demonstrates that both local and global expansion flows 
are controlled by the same physical factor - the cosmic 
vacuum with its perfectly uniform density, - when the vac- 
vnim dominates entirely. In this limit, the global cosmo- 
logical expansion and the local Hubble flow are parabolic. 
The only constant physical parameter for the both is the 
universal Hubble constant, in this asymptotic regime. 

It is also remarkable that the real Local Universe and 
the global Universe are both now rather close to this com- 
mon limit. Thanks to this underlying physical cause the 
observed values of Hi^cai and Hq are so close to each other: 
they are close because each of them is close to Hy. 

Equations 1-10 constitute the basic equations of our 
non-Priedmann model for the Local Universe. To describe 
the present-day local flow, we may put t = to — 14: Gyr 
in Eqs.1-10. Then we have the dependence of the ac- 
celeration and velocity of the flow particles on the dis- 
tance r at the present epoch. Especially informative is 
the phase (velocity-distance) relation given by Eqs.4,10, 
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showing that the velocity of a particle is determined by 
its distance and energy E, while the distance itself de- 
pends on both energy and the 'initial moment' T that are 
different for different particles (e.g. particles emerge at dif- 
ferent times from within the zero-gravity surface). This set 
of equations in combination with observational data gives 
a clear picture of the evolution and the present kinemat- 
ics of the expansion flow in the gravity/antigravity force 
field of the Local Universe. They clarify also the major 
trends in the evolution of the flow during the whole life 
time (~ 12 — 13 Gyr) of the Local group (see below). 

3. Local cosmology: General Relativity theory 

In the Newtonian theory for the local Universe (Sec.2), the 
dark energy is treated as relativistic 'fluid', in accordance 
with the General Relativity, while the ordinary matter and 
the space-time in which the matter move are described in 
terms of non-relativistic physics. Now we give a completely 
relativistic theory for the Local Universe. 

3.1. Spherically symmetric static metric 

It is most important that the set of the Newtonian spheri- 
cally symmetrical equations (Eqs.1-10 above) has an exact 
analogue in General Relativity. Indeed, there is a well- 
known GR solution for a point-like body on the vacuum 
background: 

ds^ = F{r)de - r^dn"" - F{r)-'dr'', (12) 
where 

F{r) = l-2GM/r-{r/Af. (13) 

This time-independent metric describes a spherically- 
symmetrical static space-time. Identifying here the con- 
stant parameters M and A with the mass of the Local 
Group and the vacuum timc/lcngth, respectively, we may 
use this metric for the local cosmology. The local flow par- 
ticles move along geodesies in the space-time of Eqs.12,13. 

In the limit of small deviations from the Newtonian 
(Galilean) space-time (quite appropriate for the Local 
Universe), the metric Eq.l3 takes the form: 

F{r) ~l-GM/r-^{r/A)^ = 1 + U{r), (14) 

where U{r) is the gravitational potential of the Newtonian 
theory and exactly the same as in Eq.5. 

The metric of Eqs.12,13 representing the static GR 
model of the Local Universe contains all information about 
possible motions in this static space-time. In particular, 
the radial motions of the local flow particles arc naturally 
the same (practically exactly) as in the Newtonian version 
(Eqs.1-10) of the model. 

3.2. Embedding the Local Universe into global space 

Furthermore, it is important that the 3D volume of the 
Local Universe with the metric of Eqs. 12-14 may be con- 



sidered as 'embedded' in the Friedmann global 3D cosmo- 
logical space with the non-static metric 

ds"^ = dt^ - a{t)[dn'^ + dx% (15) 

where a(t) is the 3D curvature radius and/or the expan- 
sion scale factor (for k = 1) which is proportional to the 
Friedmann radius R{t, x) of Eq.ll. 

Indeed, the Local Universe may be described as a 
spherical 'vacuole' in the uniform distribution of matter, 
as given by the exact Einstein-Straus (1945) solution. The 
size of the vacuole i?L(f, M) increases with time in accor- 
dance with the Friedmann expansion solution Rl ex. a{t), 
and the central mass M is equal to the mass of the (non- 
relativistic) matter within the sphere of radius Rl in the 
smooth matter distribution of the Friedmann Universe: 

M = ^p{t)RL{tf = Const. (16) 

In the vacuole model with the conjunction condition 

of Eq.l6, the environment of the Local Universe does not 
affect the dynamics inside it (c.f. condition 1 assumed in 
Sec.2). The Einstein-Straus (1945) solution demonstrates 
how the local non-Friedmann static cosmological model 
may be compatible with the Friedmann metrics (and dy- 
namics) of the global cosmological expansion. 

It is easy to see that the zero-gravity surface has the 
radius Rv = \{Pm / PvY'"^ RLito) ~ {^Y'^^RLito), where 
Rhito) is the vacuole radius at present. The distance in- 
terval (0.6 — l)RL{to) is dominated by the vacuum in the 
Local Universe. 

Arguing further along this line, one can imagine that 
the whole Universe may contain not one, but many (or 
inflnite number) of vacuoles of various sizes and masses. 
Moreover, a picture is theoretically possible in which vac- 
uoles fill almost all cosmic space without intersecting each 
others. Such a complex non-uniform expanding structure 
is exactly(!) described by the equations above. Obviously 
in this picture the tiny contribution of uniform matter 
distribution between vacuoles can be neglected, and one 
obtains a highly non- uniform, but completely regular, cos- 
mological model. The model describes the global expan- 
sion in terms of the relative motions of discrete masses 
which are in the centres of the vacuoles. The masses move 
apart from each others on the uniform vacuum back- 
ground, in accordance with the Hubble law and the ex- 
pansion factor given by the Friedmann theory (Eq.ll). In 
fact, because the zero-gravity radius is inside the vacuole 
in such a model, each vacuole expulses every other one 
and the expansion is generally accelerating (c.f. sect. 7. 2 
in Teerikorpi, Chernin, Baryshev 2005). 

In this idealized model, the global Hubble constant, 
H = R/R = d/a, is exactly the same for any two 
masses of the global expansion fiow. In accordance with 
the Friedmann theory (Eq.7), it is presently {t = to): 

H{to) = Ho = Hv[l + PM{to)/pvY'^ -Hv{l+^/lY'^ 
~ l.2Hv 72 km/s/Mpc, (17) 
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in good agreement (not surprisingly) with the concordance 
data (Spergel et al. 2002); it is also taken into account 
that E = in Eq.ll, according to the same data {pM = 
3 X 10~^° gcm~^ is the concordance value for the mean 
matter (dark matter + baryons) density at present). 

In terms of General Relativity, the asymptotic simi- 
larity of the local and global expansion flows is seen as 
clearly as in the Newtonian treatment of Sec. 3. Indeed, 
the static metric of Eqs.12,13 and the expanding metric 
of Eq.l5 have a common asymptotic in the limit when 
and where the vacuum antigravity dominates entirely. 
This limit is described by de Sitter's static solution. As 
is well-known, this solution has the metric of Eq.l2 with 
F{r) = 1 — (r/A)'^. The space-time of de Sitter's solution 
is determined by the vacuum alone which is ever static it- 
self, and the vacuum controls completely both space-time 
metric and any geodesic in this space-time. 

(Note that the vacuole model can be generalized: in- 
stead of the central mass one may use a spherically sym- 
metrical matter distribution of the same total mass, as 
in the Tolman-Bondi solution that is discussed e.g. by 
Gromov et al. (2001) in the context of the Local Universe.) 

4. The efficiency of vacuum cooling 

As demonstrated in Sees. 2, 3, the expansion flow described 
by the local cosmological model tends to a regular kine- 
matical structure with the linear velocity-distance rela- 
tion, f = r/A, and the universal Hubble constant Hy = 
1/A, at the limit of the cosmic vacuum domination. It 
is important to emphasize that the regular self-similar 
parabolic regime is a generic asymptotic for all the ex- 
pansion flow motions, independently of their initial condi- 
tions, in the Local Universe. 

4.1. From initial chaos to asymptotic order 

Indeed, in the model of Sec. 2, 3, the initial conditions are 
given in terms of the functions E{x) and T(x). It is obvi- 
ous from Eqs.4,9 that asymptotically both functions be- 
come negligible, since E/r^ and T/t 0, in this 
limit. The self-similar motion described by Eq.lO is a 
kind of dynamical attractor for a wide variety of motions 
possible in the gravity /antigravity potential of the Local 
Universe. If the energy condition of Eq.6 is met, any parti- 
cle escapes from the interior of the zero-gravity sphere and 
its motion reaches asymptotically the parabolic regime 
with the universal Hubble constant. This basic evolution- 
ary trend in the local flow drives it from initial chaos to 
asymptotic order (see also Chernin et al. 2004, 2005). It 
is also because this internal trend is so quick that the ve- 
locity dispersion of the local flow is rather low. 

To put this in a quantitative way, let us introduce a 
'random' velocity t; as a difference between the velocity 
f and a 'regular' (asymptotic) velocity r/Ay- Then the 
velocity v at the zero-gravity distance r = Ry is 

v{Rv) = vi = {Rv/A)[{3 + EA^Rlfl^ - 1]. (18) 
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When the trajectory of a body reaches a distance r 
outside the zero-gravity surface, the random velocity v{r) 
(for an arbitrary E) is given by the relation: 

v{r) = {Rv/A)[{r^/Rl + 2 + IRvjr + ^vxAjRy^ 

vlA^IRlfl-'-rlRy]. (19) 

Note that different values of the total energy E give an 
'initial' distribution of the velocity v. A comparison of 
v(r) with v\ enables us to see that the random velocity 
decreases along the trajectory and approaches asymptot- 
ically zero no matter what its initial value v\ is. 

4.2. The vacuum cooling factor 

With the 'vacuum cooling factor', qy = vi/v, we may give 
a measure to the efficiency of the cooling mechanism: 

qv = {viA/Rv)[{r^/Rl + 2 + 2Rv/r + 2viA/Rv+ 

vlA^R^y/^ -r/Rv]-^. (20) 

In the simplest case of the parabolic motion (£■ = 0), the 
velocity vi = (v^— l)Rv/Av, and we have: 

v{r) = {Rv/A)[{r^/Rl + 2Ry/rf'^ - r/Rv]. (21) 

Then the cooling factor in the parabolic expansion flow is 

qv{r) = v^/v{r) = {^^l-l)[{T^ IRl+2f'^-r/Rv]-\{22) 

One may see that the random velocity is diminished by 

factor 3 during the time when the body covers the path 
from r = Ry to the distance r ~ 2Rv — 3 Mpc. When the 
same body reaches, say, distances 4 or 6 Ry, the cooling 
factor increases to qv = 12 and qv = 56, correspondingly. 

Now we compare the vacuum cooling with the usual 
cosmological adiabatic cooling. The latter is described by 
the relation va(t) = constant, which gives only the factors 
2, 4 and 6 when the scale factor R{t) increases 2, 4 or 
6 times; this is considerably less than the corresponding 
factors 3, 12 and 56 in the expansion flow of our model. 

It is also interesting to follow - for a contrast - the dy- 
namics of the same expansion, but in the absence of vac- 
uum. With the same statement of the problem as above, 
but with pv = 0, we would have: 

f2 = 2GM/r + 2E. (23) 

When, for instance, .E > 0, the asymptotic (as r goes to 
infinity) motion is inertial and the velocity-distance law is 
a linear one: r = r/t, where the 'Hubble factor', H — 1/t, 
depends on time only. This means that the flow may cool 
down without the vacuum just because of its expansion. 
But the efficiency is then much lower. Indeed, taking r/t as 
a regular expansion velocity, we may introduce the random 
velocity v by the relation r = r/t + v. Then 

.= (2i^)V2[(i + ^)i/2_i]. (24) 

In the first approximation for v/(2E)^^^ < 1, this leads 
to the 'adiabatic' relation vr = constant. As a result, the 
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cooling factor in the absence of vacuum, q (x r, is signifi- 
cantly smaller than the vacuum cooling factor qy above. 

Thus, vacuum cooling is very efficient and results in a 
dynamical evolution of the outflow from chaos to order. 
When applied to the nearest distances 1-3 Mpc, the non- 
Friedmann cosmology (for parabolic expansion) gives a 
typical velocity dispersion between vi = (\/3 — l)Rv/A = 
60 km/s at r = Ry and v{r = 2Rv) = ^vi = 20 km/s 
at r = 2Rv, which is near the observed figure 30-60 km/s 
(Karachcntscv & Makarov 2001; Karachentscv et al. 2002, 
2003a; Ekholm et al. 2001; Paturel & Teerikorpi 2005). 

5. Significance of the universal Hubble constant 

In this and the next sections, we discuss new observational 
prospects offered by the local cosmology. The background 
is naturally the global cosmology where both Hy and _ffo 
are among the basic parameters. As we show here, there 
are good prospects for the uniform vacuum density and 
the universal Hy to be determined locally, but the ex- 
act value of the Hubble constant Hq refers to the average 
global expansion rate that cannot be measured in a local, 
inhomogeneous environment. Locally, a good advantage 
is the high, 10 % or better, accuracy of distance mea- 
surements, while on large scales distances continue to be 
hampered by low accuracy and systematic errors. 

5.1. The stochastic Hubble expansion rate 

Of course, a "precision cosmology" value of Hq may be ob- 
tained from the analysis of the CMBR. However, this still 
depends on assumed cosmological physics and dark mat- 
ter and dark energy components (Spergel et al. 2003), and 
it would be very important to have an independent direct 
measurement on large scales. Deep space measurements of 
Hq, even bypassing local calibrations, from the Sunuyaev- 
Zeldovich effect and the time delay in gravitational lens 
images, are still affected by rather large systematic er- 
rors. The same must be said about the proposed way of 
studying the distance scale via the Eddington efficiency of 
luminous quasars (Teerikorpi 2005). 

In intermediate and closer space, the HST Key Project 
(Freedman et al. 2001) derived i?o = 72 ± 8 km/s/Mpc, 
using secondary distance indicators calibrated from local 
galaxies whose distances are known from Cepheids. This 
value agrees with the value from the CMBR analysis and 
is a part of the concordance picture. However, one should 
not lose from sight the fact that Hq has not yet been 
determined directly on large scales with a high accuracy 
and there are still possible sources of systematic error in 
the HST value for Hq (e.g. Paturel & Teerikorpi 2005). 

Further we wish to emphasize that the measured 
Hubble constant is not the same between all pairs of galax- 
ies or clusters - contrary to the simple vacuole model. The 
expansion rate varies depending on the matter and gravity 
field variations in the space between the observer and the 
object. Therefore the really measured Hubble parameter 



is rather a stochastic variable that can be expressed rel- 
ative to the truly universal vacuum Hubble constant Hy 
and the global time-dependent average Hubble constant 
Hq = Hy{l + °^ ) • In any particular measurement 

Hi = Hv{l + EMM)y^(^i + Ai) km/s/Mpc. (25) 
Pv 

The deviation Aj must depend on the spatial scale of mea- 
surement, related to the average matter density contrast. 

5.2. The local density contrast 

The situation at small distances is also not the ideal 
one described by the vacuole model. E.g., the conjunc- 
tion condition of Eq.l6 may not be fulfilled in the Local 
Universe, unless we consider a small volume around the 
Local Group. This is because the mass of the Local Group 
is considerably less than the mass of a uniform sphere 
with the cosmological density pAiito) and the radius 5-10 
Mpc. In fact, a cosmological volume containing the non- 
relativistic mass M of the Local group (dwarf galaxies 
contribute negligibly within a few percent accuracy) has 
the present-day radius 

R{M) = [M/i^pMp^ = i?y(2-^)V3 ~ 2.2 Mpc. (26) 
o Pm 

As this is below typical distances between galaxy groups 
locally, there may be a density deficit. This would not 
influence the validity of the very local model that depends 
only on the mass of the Local Group and the value of 
the vacuum density, but it would naturally change the 
expected value of the expansion rate of the Local Universe. 

On the other hand, if we allow for an extended dis- 
tribution of dark matter beyond the 1 Mpc radius of the 
Local Group, perhaps following a fractal density law (e.g. 
Teerikorpi et al. 1998; Baryshev & Teerikorpi 2002; 2005), 
then with the fractal dimension D = 1.2 or the correla- 
tion exponent 7 = 1.8, the cosmological density would be 
reached around 3.2 Mpc. With D = 2, this distance would 
be 8 Mpc. We do not know the exact distribution of dark 
matter, but clearly it could add to the total mass so that 
the conjunction condition is locally fulfilled. 

Inspection of Fig. 4 in Karachentsev et al. (2003b) 
shows that within about 3 Mpc the Local Group is alone, 
while within 6 Mpc there are six other groups. Roughly, 
this would correspond to an average matter density of 
O.'ipM within 3 Mpc, if there is no dark matter related 
to the Local Group beyond 1 Mpc, while with the dark 
matter distribution with 7 = 1.8 the average matter den- 
sity would be close to pm- Hence, it seems that locally 
the average matter density could be, roughly, in the range 
(0.3 — 1)pmi depending on the distribution of the dark 
matter around the groups. This would agree with other ev- 
idence that within a distance of 8 Mpc there is only a small 
density excess, if any (Hudson 1992; Maccio et al.2005). 
We note that generally a randomly placed observer would 
be expected to be surrounded by a clear density excess 
(c.f. sect. 8 in Teerikorpi, Chernin, Baryshev 2005). 
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5.3. The meaning of the local Hubble expansion rate 

Because of the relatively low matter eontent around us, 
the domination of cosmic vacuum may be stronger in the 
Local Universe than on the average over the Universe as 
a whole in the present epoch. This special feature of the 
Local Universe could offer better chances for vacuum mea- 
surements there. 

First, if there is a density deficit, the local Hubble con- 
stant H in general is less than the global constant Hq, and 
closer here to the universal Hubble constant Hy, at least 
on sufficiently long scales : Hy < Hct < Hq. Secondly, 
dwarf galaxies close to the Local Group may be sensi- 
tive indicators of the vacuum, especially if the LG mass is 
concentrated within a small volume (the local cosmology 
model; sects. 2,3). We discuss this in the next section. 

To illustrate these considerations, we take a look at 
Eqs.1-10 of the local cosmology. If one assumes that the 
local flow contains mostly parabolic or nearly parabolic 
trajectories, then the very local Hubble rate is 




r/R, 



Mr) 



Fig. 1. The Hubble diagram for very nearby galaxies, with 

the units By = 1.3 Mpc and HyRv = 78 km/s have been 
used for x- and y- axis, respectively. The linear velocity- 
distance relation is shown by two straight lines, one for the 
universal Hubble constant Hy, and the other for the global 
Rv .Hubble constant Hq. The trajectories for the parabohc 

7.3 ' ' ''motion [E = 0) and for the lower limit velocity Vq are 



When 2GMAlr/r^ = 2{Ry/rf is small, the flow is 
close to the asymptotic regime, and the deviation of Hq 
from Hy is less than, say, 6% at distances r ^ 2.5i?v ~ 
3-4 Mpc. 

Of course, there arc distortions of the local flow caused 
by the Virgo cluster and other environmental effects, but 
these may still be relatively weak at r < 3 — 4 Mpc. This 
is supported by various observational studies of the lo- 
cal Hubble flow, as summarized in Teerikorpi, Chernin, 
Baryshev (2005). In these studies, the local Hubble con- 
stant is found within the interval 56 — 70 km s"^Mpc~"'^, 
so that the central value here is quite close indeed to the 
theoretically expected value of Hy. 

6. The shortest-distance indicators of dark energy 

The theory of the local flow (Sees. 2, 3) leads to two clear 
predictions that may be falsified by observations: 

(1) There is a lower limit for the flow velocities outside 
the zero-gravity sphere. Indeed, the escaped galaxies have 
the energy E > Umax = —^GM/Ry, as is seen from Eq.6. 
For any distance r, this inequality implies that the flow 
velocity must be large enough: 



r(r) > Voir) = HyRy{2Ry/r + r^/A^ - 3)^/2. 



(28) 



(2) There is a lower limit for the distances in the local 
Hubble flow: r > Ry. 

Prediction (1) is critical: if only a single isolated and 
accurately observed galaxy is found to violate this limita- 
tion, the theory or the values of its parameters must be 
rejected or modified. Prediction (2) is closely related to 
the observed kinematical structure of the flow and may 
serve for effective measurements of dark energy density. 



given. 



6.1. Comparison with observational data 

Here we concentrate on the nearest distances which are 
most sensitive to the test. It is also here where the dis- 
tances (TRGB, Ccphcids) are expected to be less in- 
fluenced by selection effects causing systematic errors 
(Teerikorpi 1997; Paturel & Teerikorpi 2005). 

Indeed, recent data on the Local group vicinity 
(1-3 Mpc from the group barycenter) are given by 
Karachentsev et al. (2002; 2004). In the critical distance 
interval 1-2 Mpc, 8 galaxies wctc studied. Their distances 
are measured with a typical reported accuracy of 10% 
and the accuracy of the velocities is 1-2 km/s. The galaxy 
SagDIG is the nearest to the Local group barycenter: its 
distance is 1.2 Mpc. Its velocity, 23 km/s relative to the 
group barycenter, is the lowest in the flow. The next clos- 
est galaxy is SexB with the distance of 1.6 Mpc and the 
velocity of 111 km/s. The three galaxies, Antlia, NGC3109 
and SexA, have almost equal distances near 1.7 Mpc and 
the velocities of 66, 110 and 94 km/s, respectively. The 
galaxy KKR25 is of the distance 1.8 Mpc and the velocity 
68 km/s. The galaxies E294-010 and KKH98 are located 
at the distance near 2 Mpc and have the velocities 81 
and 152 km/s. We also add NGC300, with its Cepheid 
distance (1.93 Mpc) recently measured by Gieren et al. 
(2004), and the velocity 125 km/s and seven other galax- 
ies with r < 2.5 Mpc from the catalogue of Karachentsev 
et al. (2004). Thus we have all the galaxies with known 
distances (from the TRGB or Cepheid method) between 
1 Mpc and 2.5 Mpc. 

All the 16 galaxies with their barycentric distances 
from the LG are plotted on the Hubble diagram in Fig.l, 
where the units Ry = 1.3 Mpc and HyRy = 78 Mpc 
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have been used for x- and y- axis, respectively. The lin- 
ear velocity-distance relation is shown by the two straight 
lines, one for the universal Hubble constant Hy ^ and the 
other for the global Hubble constant i/o- The trajectories 
for the parabolic motion {E = 0) and the lower limit ve- 
locity Vo are also given. We may see in Fig.l that all the 
galaxies here obey the restriction of Eq.28. 

We may also see that - in agreement with the discus- 
sion of Sees. 2, 3 - the two curves reveal a trend for con- 
verging to the universal (with Hy) linear velocity-distance 
relation. The lower limit velocity approaches the relation 
rather closely at the distances 5-7 Mpc (mentioned also 
in the section above) that could provide another sensitive 
interval for an observational test. 

6.2. Towards weighing dark energy 

In global cosmology, the dark energy of cosmic vacuum has 
revealed itself in observation (Riess et al. 1988, Perlmutter 
et al. 1999) at very large distances comparable with the 
vacuum length A or the horizon radius. The corresponding 
redshifts are near 



Similarly, the density of the vacuum dark energy is 



l + zv = [2pv/pM{U^)Y'^ 



(29) 



At 2; = ~ 0.7, the global cosmological expansion tran- 
sits from deceleration to acceleration. 

The rcdshift zy is a clear global counterpart of the 
local zero-gravity radius Ry'- the local flow transits from 
deceleration to acceleration at r = Ry- The difference, 
however, between them is essential. In the global cosmol- 
ogy, the balance of gravity and antigravity occurs only 
for one moment of time (7-8 Gyr ago) and simultaneously 
for the whole Universe. On the contrary, in the local cos- 
mology, the balance exists for ~ 10 Gyr, but only at the 
zero-gravity surface. In a sense, what is temporal in the 
global cosmology is spatial in the local one, and vice versa. 

Could observations at short distances near Ry — 1 — 2 
Mpc be as useful and important as the observations at 
high redshifts near zy ^ 11 

The theory prediction (2) above indicates that the dis- 
tance Ry plays a key part in the dynamics of the local 
flow. As first mentioned in Paper 1, the local Hubble flow 
starts very nearly at the zero-gravity sphere, if not exactly 
at r = Ry. This is hardly just a funny numerical coinci- 
dence with no physical sense. Rather, the coincidence is 
quite suggestive: it points out the shortest distances where 
the dark energy may be detected and measured from its 
effect on velocities. 

Indeed, as is seen from Fig.l (and from the basic equa- 
tions of the theory as well), the velocities of the flow tra- 
jectories have minima at r = Ry . The location the minima 
is the same for all the trajectories, independently of their 
initial conditions. If the location of the velocity minima 
r = Ry is found from observations, the universal Hubble 
constant may be directly determined (following from the 
definition of Eq.3): 



60 ( 



M 



I.SxIQI^Mq" Ry 



g/cm3. (31) 



Now we turn to the data of Fig.l again and inspect 
them for the velocity minimum location. Noting that the 
velocities appear to be still decreasing around r/1.3Mpc w 
1.4 a robust conservative estimate follows for the upper 
limit value of the zero-gravity radius: 



Ry ^ 1.8 Mpc. 



(32) 



This agrees with the theoretically expected provisional 
value 1.3 Mpc that we had in mind above. Combined with 
Eqs.30,31, this leads to observational restrictions on the 
values H^/M and py/M. Assuming M = 1.5 X lO^^M©, 
one has from Eq.32 an independent estimate for the lower 
limit of the dark energy density in the Local Universe: 



pyk,3x lO"'^" g/cm-^. 



(33) 



If the lowest velocity in the fiow (that of the galaxy 
SagDIG) is tentatively identified as showing the location 
of the velocity minimum in the flow, we may get Ry = 
1.2 ± 0.1 Mpc. With this (not so reliable) estimate, the 
value spread is somewhat narrower and the accuracy of 
the estimate is about 20 %: 



py = {9±2)x lO"-^" c?/cm^ 



(34) 



1.5 X 1012Mfr 



© 



R'i/ 



We give these values as examples of what could be done 
with higher precision using considerably more numerous 
distance determinations in the interval 1 - 2.5 Mpc. One 
expects many more still unmeasured dwarf galaxies within 
these distances (Karachentsev et al. 2004). In any case, 
the very initial local flow with the shortest distances and 
lowest redshifts - in drastic contrast to the important, 
but very difficult high redshift SN studies - have proved 
to be a sensitive indicator and 'weighing device' of the 
dark energy. 

7. Discussion 

The current standard cosmological model may actually be 
seen as containing two fundamental and equally important 
sectors. One is the traditional global cosmology for the 
spatial scales of 300-1000 Mpc and more, while the other 
is the local cosmology. The former studies on largest scales 
the uniform and evolving space-time; the latter deals with 
the non-uniform static space-times of the Local Universe 
and many (or infinite number) of other more or less similar 
local regions. The global sector uses the Priedmann solu- 
tion as the theory basis, while the local one needs a non- 
Friedmann theory that we proposed and described above 
(Sees. 2, 3) in a simple analytical version. We explained also 
(Sec. 3) how the two sectors of the whole model might be 
linked to each other, - at least, in our example of an ideal 
vacuole scheme. 

The omnipresent and perfectly uniform dark energy 
of the cosmic vacuum is the major physical factor which 
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unifies the two sectors of the model. It is especially impor- 
tant that, in the present-day Universe, the vacuum with 
its antigravity dominates dynamically on global scales 
and also within individual local volumes outside the zero- 
gravity surfaces. On the global scales, the vacuum anti- 
gravity accelerates the general cosmological expansion 
produced by the initial Big Bang. Locally, the vacuum 
antigravity makes dwarf galaxies ( mostly, but also galaxy 
groups outside each others' zero-gravity spheres) move 
from the Local Group and apart from each other with 
acceleration, finally establishing the Hubble fiow. 

On both global and local scales, the vacuum provides 
the motions with a common asymptotic behaviour that is 
the self-similar parabolic expansion entirely controlled by 
the vacuum antigravity. This dynamical regime is an at- 
tractor for a wide range of real motions on the cosmic vac- 
uum background, including the Friedmami uniform expan- 
sion and the local spatially non- uniform flow. In this limit, 
the expansion follows the exponential law, both globally 
and locally: r cx i? oc exp(t / Ay) ■ As a result, the lin- 
ear velocity-distance relation establishes, r cc r; R (x R, 
i.e. the Hubble expansion law. Both locally and globally, 
the expansion rate is given by the same universal phys- 
ical constant related directly to the vacuum density, in 
the limit: Hy = r/r = R/R = l/Ay ~ 60 kms-^Mpc-i. 
Correspondingly, the maximally symmetric space-time de- 
scribed by the de Sitter metric serves as the universal 
asymptotic for both Friedmann global non-static space- 
time and the non-Priedmann local static space-time. 

The concordance data for global cosmology and the 
data for the Local Universe indicate that the present-day 
states of both global and local expansion flows are not far 
from the universal asymptotic regime. It is because of this 
fact that the observed local expansion rate is so close to 
the global Hubble constant (which has been rightly consid- 
ered as a big cosmic puzzle - see Sandage 1999). Thanks to 
the same fact, the velocity dispersion is rather low in the 
local expansion flow (which has also been a true mystery). 

The simplest analytical version of the local model de- 
scribed in this paper is supported by computer simulations 
(Chernin et al. 2004, 2005) which extend the model into 
the interior of the zero-gravity surface. The simulations 
describe the Local Group as a binary (the Milky Way and 
M31) galaxy, in accordance with the Kahn-Woltjer (1959) 
classical model. They confirm that the differences from the 
point-mass analytical model are indeed small outside the 
zero-gravity surface. The binary produces relatively small 
deviations from splieric;al symmetry in the local gravity 
field at r > and also a relatively weak time depen- 
dence of the gravity field at these distances. In particular, 
the simulations show that the zero-gravity surface is al- 
most perfectly spherical at present. 

The simulations describe all the history of the local 
flow during 12-13 Gyr since the formation of the Local 
Group. They enable us also to study the "initial condi- 
tions" for the local flow, and indicate that the Little Bang 
model (Byrd et al.l994) for the initial chaotic state of the 
Local Group (see also van den Bergh 2003) could produce 



a flow of galaxies across the zero-gravity surface and their 
radial velocities are later amplified by the antigravity push 
of the cosmic vacuum outside the zero-gravity surface. 

Our analytical theory and simulations are in good 
agreement with the large N-body ACDM cosmological 
simulations (Ostriker and Suto 1990; Suto et al. 1992; 
Governato et al. 1997; Nagamine et al. 1999, 2000; Maccio 
et al. 2005). These simulations give important insights 
into the early history of the Local Universe when its mat- 
ter started to separate from the general cosmological ex- 
pansion at the cosmic age of 0.5-1.5 Gyr. The picture is 
complementary to ours and demonstrates that the Local 
Universe with a massive group in its center and the quies- 
cent expansion outflow is rather typical for scales of a few 
megaparsecs. Such local structures emerge in a natural 
way from the standard picture of gravitationally insta- 
bility with the Harrison- Zeldovich initial spectrum. Very 
importantly, these 'global' simulations have shown con- 
firming the original approach of our Paper I that the 
dark energy is necessary for explaining cool local fiows. 

8. Conclusions 

We summarize briefly the main results of the present 
study: 

• A non-Friedmann model with a non-uniform static 
space-time was proposed for the neighborhood of the 
Local Group in order to describe the trajectories of 
galaxies flowing across the local zero-gravity surface 
at r = Ry — 1-3 Mpc. The model is given both within 
Newtonian theory and General Relativity. 

• It was shown that this local space-time, when treated 
within General Relativity, may be embedded in the 
global Friedmann space-time, and the Local Universe 
may described as a spherical vacuole, according to the 
Einstein-Straus solution. 

• We emphasize the difference between the truly univer- 
sal Hubble constant Hy, depending just on the con- 
stant vacuum density, and the measured local or global 
Hubble constant Hq that depend both on spatial mat- 
ter distribution and scale and on the cosmic time. 

• We have inspected the predicted trajectories of the 
out-flowing galaxies and demonstrated analytically the 
efficiency of vacuum cooling, which is higher than the 
usual adiabatic cooling: all trajectories within the vac- 
uole tend towards the Hubble relation with the univer- 
sal vacuum Hubble constant Hy. This gives a natural 
explanation for both the low velocity dispersion in the 
flow and its closeness to the global Hubble constant. 

• Still another problem, the absence of turning-down of 
the velocity-distance relation close to the Local Group 
(Ekholm et al. 2001; Thim et al. 2003) and in general 
the absence of the infall towards local groups (Whiting 
2005) may be explained by the present model, where 
the out-flowing dwarfs remain close to the universal 
Hubble relation and have a positive minimum velocity 
at r = Ry. 
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• An interesting result is that all out-flowing trajectories 
have a velocity lower limit at any distance, correspond- 
ing to the total energy just required to reach r = Ry 

• Preliminary comparison with available distance and 
velocity data for nearby galaxies, 1-2 Mpc, shows 
agreement of the predicted trajectories with E < 
and may suggest the expected minimum around 1.3 
Mpc. The local Hubble diagram V/{HvRv) vs. r/Ry 
is a promising instrument for such studies. 

Finally, the local Hubble flow is not so mysterious any 

more, but it may rather be seen now as a natural tool 
for the detection and study of the vacuum antigravity and 
the dark energy. In fact, we have already used it for these 
purposes since our Paper I. Combining our theory and 
computer simulations with precision observations, we have 
discovered dark energy in the Local Universe. 

• The flow provides strong evidence for the very fact of 
the presence of dark energy in the Local Universe on 
the small scales down to ~ 1 Mpc; the local DE den- 
sity seems to be close to that measured on the largest 
cosmological scales. 

• The dark energy in the Local Universe can be well 
described in the simplest way as the cosmological con- 
stant or cosmic vacuum with the dcnsity-to-prcssure 
ratio equal to -1; this agrees with the concordance 
global cosmological model. 

• The dark energy with these properties might have ex- 
isted unchanged for 12-13 Gyr in the Local Universe, 
since the formation of the Local Group; the global cos- 
mology has no data on the dark energy at such dis- 
tances/times. 

These conclusions need and deserve further observa- 
tional verification. We advocate for extensive systematic 
high precision observations of galaxies, focused on short 
distances, especially of 0.5-2 Mpc. The advantages of such 
data are clear. These observations may give us increas- 
ingly exact information about the dark energy. In partic- 
ular, the local measurement of the dark energy density 
can be much improved as compared with our preliminary 
estimate here. Strong constraints may also be obtained on 
(1) deviations from the standard density-to-pressure ratio 
in the present-day Local Universe, (2) the time variation 
of the dark matter density for a period of 12-13 Gyr since 
the formation of the Local Group, and (3) the time varia- 
tion of the density-to-pressure ratio during the same long 
cosmological time. 

Observations near the zero-gravity surface of the ra- 
dius Rv ^ 1 Mpc may prove to be as useful as the observa- 
tions near the redshift 2y ~ 1 of the zero-gravity moment 
in the global Universe. With increasingly precise measure- 
ments of distances and velocities in the Local Universe, the 
local sector will be a central arena of further progress in 
cosmology. 
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